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^ . Abstract 

CN , We explore the addition of fundamental matter to the Klebanov-Witten 

^ \ field theory. We add probe D7-branes to the M = 1 theory obtained from 

\ placing D3-branes at the tip of the conifold and compute the meson spectrum 

I for the scalar mesons. In the UV limit of massless quarks we find the exact 

dimensions of the associated operators, which exhibit a simple scaling in the 
large-charge limit. For the case of massive quarks we compute the spectrum of 
scalar mesons numerically. 



1 Introduction 

The gauge theory/string theory correspondence [H [21 [3] furnishes a powerful set of 
tools for understanding gauge theories at strong coupling by performing computations 
in a dual string theory at weak coupling. However, the correspondence is only well- 
understood in systems where the string background is highly symmetric and nearly 
flat, but we expect that the duals to many interesting gauge theories (such as large- 
QCD or SQCD) will not have these properties. It is therefore an interesting challenge 
to study less symmetric string backgrounds, and in particular to study backgrounds 
with reduced supersymmetry. 
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One interesting class of models arises from compactifications of string theory on 
noncompact Calabi-Yau manifolds with D3-branes at conical singularities, which 
generically give rise to A/" = 1 gauge theories with product gauge groups and bi- 
fundamental matter. These models are attractive for several reasons. They possess 
minimal supersymmetry and are therefore closer to realistic gauge theories than the 
well studied = 4 case; also, they lead to conformal field theories where the quantum 
conformal invariance is not obvious by inspection of the field theory (but where the 
supergravity dual makes conformal invariance manifest.) Perhaps the most striking 
feature of these theories is that one can break conformal invariance in a controlled 
way by adding fluxes through cycles of the Calabi-Yau geometry, which induce RG 
flow and confinement at low energies. 

However, one missing element of these models is fundamental matter. Aside from 
being experimentally important, fundamentals give rise to many interesting things 
such as the phase structure of super- Yang- Mills theory in the infrared. In confining 
theories the fundamentals of course do not appear as asymptotic states but are instead 
confined in mesons and baryons. 

In this note we study the mesonic fluctuations of a particular set of mesons in 
the conifold theory of Klebanov and Witten P|. This theory is interesting for its 
relative simplicity and also because its non-conformal version flows to a theory very 
similar to A/" = 1 pure glue theory in the infrared. Moreover all metrics for the 
corresponding supergravity solutions are known, allowing explicit computations. The 
mesons which we study arise as fluctuations on D7-branes which are embedded in 
the string background. The fundamental fields come from strings connecting the 
stack of D3-branes to the D7-branes. In the usual decoupling limit, the 3-7 strings 
and 3-3 strings, which describe the gauge theory, have a dual description in terms of 
the closed strings and 7-7 strings. The closed strings are the usual glueballs of the 
strongly coupled field theory while the open 7-7 strings are naturally identified with 
the mesons. 

We will compute the spectrum of operator dimensions, which, as we will see, can 
be done exactly for a large portion of the states, and we will study the effect of giving 
masses to the quarks (which requires numerical work). 

The paper is organized as follows. In section [2] we review the geometry of the 
conifold. In section [3] we discuss adding flavor to the Klebanov- Witten field theory 
by the addition of probe D7-branes. In section H] we compute the spectrum for scalar 
mesons. In the case of massive quarks, we compute the mass spectrum numerically, 
but in the massless case (corresponding to the UV limit of the gauge theory) we 
obtain the spectrum analytically. In section [S] we discuss our results. 
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2 Review of the Conifold 



In this section we briefly review the geometry of the conifold in order to fix notation. 
Useful references are H El El E El E] . 

The conifold is a non-compact Calabi-Yau 3-fold, defined by the equation 

Z1Z2 - Z2,Zi = (1) 

in C^. Because Eqn.([T]) is invariant under an overall real rescaling of the coordinates, 
this space is a cone, whose base is the Einstein space T^'^ [H The metric on the 
conifold may be cast in the form [S] 

dsl = dr^ + r^ds^Li , (2) 

where 

ds^i,i = - + ^ cos 9id(f),^ 6 ^ ^^^'^ sm^9id(j)^) (3) 

i=l i=l 



is the metric on T^'^. Here ip is an angular coordinate which ranges from to An, 
while (6*1, (pi) and {62, ^2) parametrize two S^s in the standard way. This form of the 
metric shows that T^'^ is a U{1) bundle over x S^. 

These angular coordinates are related to the Zi variables by 
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It is also sometimes helpful to consider a set of "homogeneous" coordinates Aa, Bb 
where a, 6 = 1, 2, in terms of which the Zi are 

Zl = AiBi, Z2 = A2B2, (5) 
Z3 = A1B2, z^ = A2Bi. (6) 

With this parameterization the Zi obviously solve the defining equation of the conifold. 

We may also parameterize the conifold in terms of an alternative set of complex 
variables Wi, given by 

Zi=Wi+ iW2, Z2 = Wi- iW2, ^^.^ 

Z3 = -W3 + iw4, Z4 = W3 + iwi. 
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The conifold equation may now be written as 

and we identify the T^'^ base of the cone as the intersection of the conifold with the 
surface 

J2\w^\' = r\ (9) 

T^'^ described in this way is exphcitly invariant under 5*0 (4) ^ SU{2) x SU{2) rota- 
tions of the Wi coordinates and under an overall phase rotation. Thus the symmetry 
group of T^'i is SU{2) x SU{2) x f/(l). 

An important fact about T^'^ is that it has Betti numbers 62,^3 = 1- The corre- 
sponding two-cycle and three-cycle may be expressed in terms of harmonic differential 
forms: 

UJ2 = ^(^^11-^^22), (10) 

Us = CA^2- (11) 

In this paper we will consider D7-branes in the model of Klebanov and Witten P]. 
This model is a particularly simple Af = I gauge/gravity dual, obtained by placing a 
stack of N D3-branes near a conifold singularity. The branes source the RR 5-form 



flux and warp the geometry: 

dslo = Kr)-^'Mx^dx^ + h{rf'\dr^ + rMsl,,,) (12) 

h{r) = 1 + ^ (13) 

QsF^ = d^x Adh-^ + *{d^x Adh-^) (14) 
97 

= flngsNa'^ (15) 



Hereafter, we specialize to the near- horizon limit r/L <^ 1, and set L = 1 for conve- 
nience. It may be easily restored by dimensional analysis at any point. 

The dual field theory has gauge group SU (N) x SU{N) and matter fields Ai^2, -Si, 2 
which transform in the bifundamental color representations (N, N)^ and (N, N)^. The 
theory also has a superpotential 

W = \Tr{A,B,AkBi)e'h=K (16) 

By solving the F-term equations for this superpotential, we obtain supersymmetric 
vacua for arbitrary diagonal 2 and -81,2, so that the moduli space of the field theory 
is precisely that of N D3-branes placed on the conifold. 
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3 Adding flavor 



In this section we review the procedure of adding flavor branes to AdS/CFT in general 
and make several useful comments on adding flavor to the Klebanov-Witten field 
theory both in terms of the bulk geometry and the dual field theory. This general 
procedure was first pointed out in [lOl [11], [12] and was exploited in the AdS^ x 5*^ 
case in [13]. Some other examples of flavored theories with probe branes have been 
studied in[Tlll5l[l6l[I3[T8l[l9l[20l[ai[22l[23l[2l|2^ 

One way to add flavor to AdS/CFT is to take a system of D3-branes and then 
to add D7-branes which fill the four directions and four of the six transverse 
dimensions [TT]. In flat space such a configuration of branes is clearly super symmetric. 
As usual there is an A/" = 4 SU{N) SYM theory living on the D3-branes. Strings 
with one end on a D3-brane and one end on a D7-brane couple to the fields of the 
D3-brane gauge theory as quarks. 

For AdS / CFT purposes we can now take the supergravity approximation in which 
D3-branes are replaced by an AdS^ x geometry with Ramond-Ramond flux, while 
we retain the D7-branes as probes which fill the five AdS directions and which wrap a 
topologically trivial 3-cycle of the internal 5-manifold (for example an submanifold 
of the of AdS^ x S^). The triviality of the 3-cycle guarantees that the brane carries 
no net charge and will not introduce any tadpoles. On the other hand, topological 
triviality also suggests that one might be able to shrink the 5*^ and slip it off of the S^, 
naively in contradiction with the flat space picture of D3 and D7-branes. It turns out 
that subtleties of the AdS geometry play a key role in ensuring stability. The mass 
eigenvalues of modes controlling the D-brane slipping off the 3-cycle are negative, but 
are above the Breitenlohner-Freedman bound [27], so that the 7-brane embedding is 
stable. 

In the flat space picture, if the D3-branes and D7-branes intersect then the quarks 
are massless, and if the D3-branes and D7-branes are separated then the quarks are 
massive. This translates nicely into the AdS picture in the following way. A D7-brane 
which intersects the D3-branes in flat space gets mapped to a D7-brane which fills the 
whole AdS space and wraps a three-sphere of constant size in the 5*^. On the other 
hand, a D7-brane separated from the stack of D3-branes maps to a D7-brane which 
wraps an with some asymptotic size at large AdS radius, but this shrinks to 
zero size at some finite radius (which is possible because of the topological triviality). 
In the 5-dimensional AdS space the D 7-brane appears to fill out the radial direction 
up to some minimal radius where it "ends." 

It is interesting of course to consider theories with branes in spaces which are not 
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flat. The basic picture of D3 and D7 branes contributing gauge fields and quarks will 
not change, but many details are different. For simplicity throughout this paper we 
specialize to the case of a single D7-brane. If the number of D3-branes is large then 
the D7 backreaction can be systematically neglected and it is appropriate to treat the 
D7-brane as a probe, which we do throughout this paper. Inclusion of backreaction 
effects in other geometries has been explored in (TUl ESI EH] ■ 

Let us consider D7-branes embedded in the geometry of the conifold by the equa- 
tion Zi = fi. In terms of the standard coordinate system, 

3/2 W-<^i-<^,) 

2 2 

so the embedding equation gives two conditions, one on the magnitude of zi and one 
on the phase: 

2/3 

-0 = I ■ :rL 0. 1 . (17) 




ipo = 4>i + 4>2 + const. 



This embedding can be explicitly shown to be supersymmetric by considering the 
K-symmetry on the worldvolume of the brane [30]. A slightly different embedding 
equation was studied in the warped deformed conifold by [31]. 

It was proposed in [32| that the embedding zi = fi leads to fields, summarized in 
Table [1] and a superpotential of the form 

W = Wfi 

avors ~l~ ^^massesi (1^) 
W flavors = hqAiQ + gqBiQ, Wmasses = + ^2QQ- (20) 

To relate this superpotential to the D7-brane geometry, let us probe the space with 
a single D3-brane, which corresponds to giving some expectation values to Ai and 
Bi. One then finds that the theory on this probe has a massless flavor when AiBi = 
fJ^ifJ^2/{gh), which is exactly of the form of the embedding equation zi = fi. Part 
of the motivation for this superpotential was a comparison with a type IIA brane 
construction [33] where a D6-brane splits on an NS5-brane, contributing two flavor 
branes and correspondingly two sets of flavors. For the type IIB picture, in the mass- 
less limit of the fleld theory, this corresponds nicely to the presence of two solution 
branches of Zi = 0, namely = and 62 = 0. If the quarks are massless there is an 
SU{K) X SU{K) flavor symmetry, where K is the number of probe D7-branes. If the 
quarks are massive then the two branches of the D7-branes connect and the flavor 
symmetry is broken down to the diagonal SU{K). 
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Field 


T T / 7\ T \ T T ( 7\ T \ 

SU (iV,) X S\J (iVj 


SU(K) X SUiK) 




(N,l) 


(K,l) 


Q 


(N,l) 


(1,K) 


Q 


(1,N) 


(1,K) 


Q 


(1,N) 


(K,l) 



Table 1: Representation structure of the added M —\ flavors. 

An alternative perspective is to suppose that one of the masses //j is larger than 
the other and then integrate out the associated flavors. Then one obtains a quartic 
superpotential of the form 

W = q{A^B^-lx)q (21) 

which again produces the appropriate massless locus for a D3-brane probe. Our probe 
calculations will show that this quartic superpotential is consistent with adding D7- 
branes with massive flavors (and then with the limit where we take the masses to 
zero.) Of course, because wc believe the quarks can be massive the consistency was 
virtually guaranteed. However, taking the limit /i — > and setting /i = are different 
things, and it is unclear whether the theory corresponding to the cubic superpotential 
can be realized or not. 



4 Scalar mesons 

In this section, we compute the dimension and mass spectra of the scalar mesons. As 
discussed in the introduction, in the probe and decoupling limits the 7-7 strings arc 
identifled with the mesons in the dual fleld theory. We will thus be able to extract 
the mass spectrum of the spin=0 mesons and their conformal dimension in the UV 
limit by studying the 7-7 strings. 

The semiclassical dynamics of this D7-brane are captured by the Dirac-Born-Infeld 
action 

Sdbi^t, J d'^^ - det {qmn + Fmn) + ^ / C4AF2AF2. (22) 

where ^'^ are coordinates on the D7-brane. We will compute the spectrum of fluctua- 
tions for the D7-branes using this action. 
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Let us consider the fluctuations of scalar modes alone, with all D7-brane gauge 
fields turned off. Then the DBI action is simply the worldvolume of the 7-brane. Let 
us choose as coordinates on the brane eight of the spacetime coordinates: (x'', 9i, 62, 0i, 02)- 
The fluctuations can be described by setting 



The unperturbed induced metric takes the form 



(23) 
(24) 
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cot' 



\ cot f cot f 



g cot -j- cot 2 



9<i>i<i>j 

1-2 ^ 1 

\-\-\ cot^ Y 

\ sin^ ^1 + i(l + cos Oif i(l + cos ^1) (1 + cos 62) 
^ |(l + cosei)(l + cose2) Isin'^s + Kl + cos^a)' 

One expands about this metric via the matrix identity 



(25) 

(26) 
(27) 



Vdet A + M = VditA(l + ]p^A~^5A+^{TYA-^5Af - ^A-^5AA-^5A+ ...{2%) 

The terms first order in the fluctuations x and r] turn out to be total derivatives, 
as is necessary for our embedding to be a solution of the equations of motion. The 
quadratic order fluctuations lead to an action of the form 



= T7 y d^xd6id92d(f)id(t)2'\/ del qq— 

-A-0^xd<p2V - ^ ^2 02 ( ^^^1 + ) xd4>2V (29) 

sm o sm \ / 



+ 



+ 



sm 



2 01 



with 



C = 1 + - cot^ — + - cot^ — . 
3 2 3 2 



(30) 



4.1 The UV/massless limit 

Even though the quarks have mass, if we flow to the UV, they effectively become 
massless and conformal symmetry is restored, at least at the classical level. It is 
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interesting to inquire what the dimensions of the operators are in the UV field theory. 
In the dual, this corresponds to computing near the boundary of AdS. 

Examining f|T7j) we see that there are two ways to go near the boundary: ^ or 
62 ^ 0. Which one we choose will determine which side of the conifold we are on near 
the boundary. In the current setup, the physics is symmetric between exchange of 61 
and 62 so we will simply choose the ^ limit. We can compute the dimensions of 
the operators in the field theory by examining the scaling of the 7-7 strings near the 
boundary. 

We define tq = /i^/^e"^/'^ (now a good coordinate because of the conformal invari- 
ance), and cos^ y = x. Defining the linear combination of fields 

^^ = X±iv (31) 

we find that the equations for are two fully decoupled partial differential equations. 
This equation is solved by a separation of variables ansatz 

$± = p±(x)e'='^/=^e*™^'^^+^'^^'^^ (32) 

The equations of motion for the scalar reduces to ordinary differential equations for 
the functions p{x)^ which take the form 

. d d 1,,, , m?3 — x rriimo 
{l-x)—x— + -k{k~l)~^- + 



dx dx Q 8 1 — X 2(1 — x) 

ml mil+x m2 \ ±, . ^ .^ox 
=F-r^^ ±7T7^ 7 P (a;) = 0. (33) 



4x(l-x) 4 1-x 2(1 -x)^ 

This equation has singularities only at x = 0, 1, 00 and is therefore of hypergeometric 
type. To see this explicitly we define new functions by rescaling the by factors of 
X and (1 — x), which allow us to remove the terms in ([3SD proportional to ^ and 
Explicitly, we write 

p±(a;) = xP{l-xff^{x) (34) 
for which the equation of motion becomes 

x(l - x)—-^ + [l + 2p- x(l + 2p + 2q)) — + 



dx'^ dx 1 — X 



+- — 7^ + 6^^^"^)" 4 (2 ^""^ ]-ip + q) 



/±(x) = 0.(35) 



Appropriate choice of the parameters p and q eliminates the terms proportional to 
1/x and 1/(1 — x). The wavefunctions are given by 

f^{x) =2^1 (-a, a + 2{p + q),l + 2p; x) (36) 
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which are regular when a is a non-negative integer; it turns out that the original p 
are also regular with this condition over the range 0<a;<l(0<^2<7'") which 
encompasses our domain. We also find that there are two possible values of k: 

ki = ^ + ^ ^1 + Sm? ^ 6mi + 2A{a +p + qf (37) 
k2 = ^ - ^ ^1 + Sm? ^ 6mi + 24(a +p + qy. (38) 

To be painfully exphcit, we exhibit the solutions for /+ (the /~ are straightfor- 
wardly related.) It is clear that there are always two choices of p and q which do the 
trick; to make regularity transparent we will always choose p and q to be positive. 
We then have four cases: 

• 1712 > 0, nil > We choose p — 7712/2 and q = 1 + |(mi — 1712), so that 
/"^ =2-^i(— Q;,2 + mi + a,l + 777.2; x). The wavefunction is regular if a is a 
non-negative integer (negative a gives irregular or redundant solutions) and so 
the two values of k are quantized to be 



k ^ ^ ±\Jl + 3ml - 6mi + 24{a + ^ + ly. (39) 

1TI2 ^ 0,mi < 7712: Again we choose p = 7772/2, but now to make regularity 
obvious we take q = (7772 — mi)/2, such that g > 0. Now f'^ =2-Fi(— a,2m2 — 
777i + a, 1 -|- 777.2; x), again with a a non-negative integer. The quantized values 
of k are 



k = ^±^^1 + 3ml - 6mi + 2A(a + 7772 - ^)'. (40) 

7772 < 0, 777i > 7772: Now wc choosc p — —7772/2 and q — 1 + |(777i — 7772), finding 
that /"•" =2-^1 (—a, 2 + 777i — 27772 + a, 1 — 7772; x), with a a non-negative integer. 
The allowed values of k are 



k = ^ ±^J 1 + 3ml - 6mi + 24(a + ^ - 7^2 + 1)=^- (41) 

777.2 < 0,7771 < ^2- Now we choose p = —7772/2 and q = (7772 — mi)/2, finding 
that —2Fi{—a, a— mi, 1— 7772; x), with a a non-negative integer. The allowed 
values of k are 



A; = ^±^W1 + 3777? - 67771 + 24(q;-^)2. (42) 
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To find the dimensions of the operators, we recall that e^^^^ ~ r^^. In the 
AdS/CFT correspondence a minimal massless scalar field dual to an operator of di- 
mension A scales as for its normalizable part and r^~'^ for its non-normalizable 
part. However, by examining fl29l) we see that the kinetic terms for these scalars are 
not canonically normalized, which means that the possible scalings at infinity are 
modified to r~^+^ and r^"^"^^ for some p. Using the values for ki and ^2 we have 
2 A — 4 = fci — A;2, which one can compute straightforwardly. 

The dimensions are mostly complicated irrational numbers (reminiscent of the 
closed string spectrum on T^'^ [HI [31|) but a few features of the spectrum stand out. 
The lowest mode has mi = m2 = a = 0, and is simply a constant; it can be assigned 
dimension 5/2 or 3/2. From the earlier discussion of the massive field theory, it is 
natural to choose dimension 3/2 and associate this mode with the operator qq. Note 
also that the mode with mi = m2 = 1 and a = has dimension 3, appropriate for a 
superpotential term - we identify this mode with the operator qA2B2q. If added to 
the superpotential, this operator would change the D7-brane embedding from zi = fi 
to zi + ez2 = fi. 

For large mi, all the dimensions scale as A ~ 3mi/2. This is consistent with 
identification of the corresponding gauge theory operators as 

q{AB){AB)...{AB)q (43) 

where, ignoring the q fields, each insertion of (AB) should increase the dimension by 
3/2 and the relevant SU{2) charge (associated with mi) by one unit. Unlike the case 
of baryonic operators on the conifold, where one finds an exact scaling A ~ SN/A 
|35] . we see that the mesons only exhibit a simple scaling with the charge in a large- 
charge limit. For small charges there are boundary effects due to the quarks which, 
at least in the large- iV limit, are completely calculable here. This behavior should 
also be contrasted with the case of flavors added to A/" = 4 super- YM theory, where 
the meson dimensions were pure integers. 

If we take instead the limit of large a, we see that the dimensions scale as A ~ 
\/6a. It would be interesting to flnd an explanation for this curious scaling in the 
fleld theory. 

4.2 The mass spectra 

Having obtained the conformal dimensions for the mesons in the conformal limit, we 
would like to compute their spectrum by solving the full differential equation without 
taking any simplifying limits. Unfortunately, we will find that the equation is not 
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amenable to analytic solution and so we will have to appeal to numerical methods. 
We will display selected results from several cases that are illustrative of the general 
behavior. 

We will again find that the linear combination of fields flHTl) decouples the equations 
of motion. Rewriting the action ( l29l) in terms of $^ and varying gives 

^ da(^V^o9fdt^^)±3t\^=doXV^ol'''^^)~l^^]d^,^^ = 0, (44) 



where the a, h indices run over the and the 6'i^2, 0i,2, and 

7*' ^ (46, 

36 sm fc'j 

The inverse components of the metric are straightforward to find from the form given 
in ( l25ll271) . Since dx^- and d^. are Killing vectors we can write 

We find that (HH) becomes 



Note that for massive modes k"^ = k^k^ < since it must be timelike. Using simple 
Kaluza-Klein arguments we see that the mass of the mesons in the dual field theory 
is = —k"^. It is evident from the form of this equation that the only difference 
between the equation for and for ip~ is in a term proportional to rrij. Thus, we 
can choose to solve for without loss of generality. This equation cannot be solved 
analytically. In addition, we have found no simple way to separate the equation in 
the 6^1, 6^2 directions and so we must use a numerical approach to solving the partial 
differential equation. 



4.2.1 The numerical approach 

Because we are unable to separate the partial differential equations into ordinary 
differential equations, we must use a technique slightly more involved than the regular 
finite difference scheme and shooting technique. We will make use of the finite element 
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Figure 1: Mass eigenvalues versus eigenvalue number for the zero node modes with 
mi = m2 = 0, = .02. +s denote actual eigenvalues and the solid line is a best 

fit line with equation -k^ = 240^^ - Q8n + 290. 

method and the Arnoldi algorithm via Matlab to solve for the mass eigenvalues, —k"^. 
We will use a mesh with 2779 nodes and 5392 triangles for all problems involved. The 
^jS are part of two different S'^s and thus range from to vr. Since we already know 
9i ^ corresponds to going near the boundary, and we want normalizable modes, 
we place Dirichlet boundary conditions at 9i^2 = 0. We must also demand regularity 
at di 2 = IT, which corresponds to placing Neumann boundary conditions at 612 = tt. 

We will first examine the simplest case, when mi = m2 = 0. Setting = .02 

we solve (jH]) for the first 50 eigenvalues. The eigenvalues break up into different 
series corresponding to the number of nodes in the (^1,^2) plane. In figures [1] and [2] 
we display the first two such series. Higher series have similar behavior. The + signs 
denote actual mass eigenvalues, while the solid lines are best fit lines. 

We will also find similar behavior for modes with mi^2 7^ 0. In figure [3] we display 
the zero node modes for the case mi = 1, m2 = 2 with /i"^/'^ = 2 (note: changing /i 
just changes the eigenvalues by an overall scaling, as expected since it merely scales 
the mass gap for the quarks). We find similar behavior for other values of mi^2 and 
different number of nodes. 

For all the cases we see that in the large n limit we have M ~ n^i'^^^ as we would 

2 /3 

expect. Restoring L by dimensional analysis, and using nig ~ we find the the 
mass gap for the lightest meson is 



M, 



rriq 



(49) 



gap 
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Figure 2: Mass eigenvalues versus eigenvalue number for the one node modes with 
mi = 77^2 = 0, = .02. +s denote actual eigenvalues and the solid line is a best 

fit line with equation -k'^ = 27071^ + 490n + 1500. 
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2 4 6 8 10 12 



Figure 3: Mass eigenvalues versus eigenvalue number for the zero node modes with 
mi = 1, 1712 = 2, /i""^/^ = 2. +s denote actual eigenvalues and the solid line is a best 
fit line with equation -A;^ = 3.2n^ + 5.4n + 29. 
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Therefore, in the supergravity regime where ggN ^ 1 we find the meson mass is 
much smaller than the quark mass. At large t'Hooft coupling we find that the binding 
energy of the mesons almost completely cancels the rest energy of the quarks. This 
is similar to the situation in AdS^ x 5*^ [13j . 

5 Discussion 

In this note we have computed the spectrum of mesons in an A/" = 1 field theory 
corresponding to fiuctuations in the position of a holomorphically embedded D7- 
brane. In the limit of nearly massless quarks, the field theory is classically conformal, 
and also conformal at large- A^, and the spectrum turns out to be computable exactly, 
where the dimensions in general are complicated irrationals. 

There are a few operators for which the exact results are simple. Among these 
are the lowest mode, corresponding to a mass term, with dimension 3/2, and a mode 
corresponding to a BPS fiuctuation of the D7-brane, with dimension 3. The existence 
of these operators suggests that a consistent superpotential for our fiavored theory is 



It would be interesting to study the Klebanov-Strassler theory [36j obtained at the 

end of the duality cascade with the addition of 3-form fiux with this superpotential. 

In the strictly massless limit, = 0, it is possible to relax our embedding con- 
dition slightly. Specifically, with a nonzero mass we imposed a relation between the 
azimuthal coordinates, ~ 0i ~ 02 = 0. However, when the mass is zero this con- 
dition need not apply; it would be nice to see what relaxing this condition would 
mean for the field theory (in particular, whether it is possible to realize the cubic 
superpotential discussed in section 3.) 

The appearance of irrational dimensions is not surprising, in light of similar results 
for the glueball spectrum of the conifold [H [31] • However, this feature of the meson 
spectrum differs from the A/" = 4 case, where the meson dimensions were pure integers. 
In particular, we do not find a tower of states with spacing 3/2, except in the large R- 
charge limit; more precisely, in this limit the spectrum is of the form 3k/2+0{l/ J). It 
might be possible to compute these 1/J corrections in a plane- wave limit, or perhaps 
in some other formalism. It would be interesting if such a comparison with our exact 
results were possible. 

We have also numerically computed the spectrum for the case of massive quarks. 
In the large ggN limit the meson mass gap is significantly smaller than the quark 



W = qAiBiq. 



(50) 




15 



masses. We have uncovered a relatively simple quadratic scaling behavior for the 
meson masses. It would be nice to find, either with analytical or more numerical 
work, the exact functional dependence on n, mi, 1712 etc. 

All of our calculations have been in the probe limit and further studies of the 
backreaction would be interesting, especially for the Klebanov-Strassler deformed 
conifold theory. However, it may still be possible to learn things from further study 
of probe theories. In particular, it would be interesting to study the dynamics of 
nontrivial classical field configurations in the D7-brane worldvolume. Such fields 
would correspond to dissolved D3-branes or anti-D3-branes. The anti-brane case is 
particularly interesting, as it would break supersymmetry along the lines of the KKLT 
scenario j37]P1 . but with the possibility for some moduli to be fixed by the D7-brane. 
We leave these suggestions for the future. 
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